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The spectral properties of a disordered electronic system at the metal-insulator transition point 
are investigated numerically. A recently derived relation between the anomalous diffusion exponent 
n and the spectral compressibility \ aT - the mobility edge, \ ~ fl/^d, is confirmed for the integer 
quantum Hall derealization transition. Our calculations are performed within the framework of an 
unitary network-model and represent a new method to investigate spectral properties of disordered 
systems. 
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Since the development of the classical Random Matrix 
Theory (RMT) by Wigner, Dyson, Mehta and others the 
statistics of energy levels of complex quantum systems 
have become the subject of research in many areas of 
physics, from the study of atomic nuclei to the investi- 
gation of disordered metals and quantum chaos Q . The 
remarkable property the level statistics of these systems 
have in common is their universality, regardless of the 
microscopic details of the particular system. 

In the following we will focus on the spectral proper- 
ties of systems undergoing a metal-insulator transition. 
In this case at least two entirely different types of level 
statistics are involved. In the insulating phase the proba- 
bility density of an eigenstate is almost completely local- 
ized within a comparatively small volume. As long as the 
spatial extension of any two states is small compared to 
their distance, they are independent in the same way as 
two states of two separate systems. Consequently, in the 
localized regime the energy levels are uncorrelated and 
therefore governed by Poisson statistics. For example, 
the probability to find two consecutive levels separated 
by an energy e = sA, A being the average level spac- 
ing, is given by the Poisson-distribution P(s) — exp(— ,s). 
The number variance £2^) = ((n — {n)) 2 } of an en- 
ergy interval which on average contains N = (n) levels is 
£2 (N) = N, according to the central limit theorem. 

In the metallic or dclocalizcd phase, on the other hand, 
the eigenstates are extended over the entire system. In 
this case the disorder potential causes the levels to repel 
each other. This level repulsion leads to a considerable 
rigidity of the spectrum with respect to fluctuations in 
the level density: the number variance increases only log- 
arithmically with the number of levels, £2 (AT) oc ln(A) 
and the spectral compressibility defined by 
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vanishes; in the delocalized regime the spectra are in- 
compressible. The level spacing distribution is described 
by the so called Wigner surmise, P(s) oc exp(— C/3S 2 ), 
where is a number of order one which depends on the 
symmetry of the system's Hamiltonian Q. The factor 
s 13 , missing in the Poisson statistics, reflects the strong 
level-repulsion. 

At the mobility edge, where the two phases with their 
different kinds of spectra meet, things become more com- 
plicated and have given rise to extended investigations 
and controversial discussions concerning the critical level 
statistics . In the vicinity of the transition energy 
the probability density of a state is neither localized on 
a small confined area nor smeared out almost homoge- 
neously over the whole system, but forms a self-similar 
measure which fluctuates very strongly on all length 
scales. It is best described in terms of multifractality 
Q. Since the repulsion between levels is strongly influ- 
enced by the spatial correlations of the corresponding 
eigenstates, their multifractal structure enters the level 
distribution. 

In recent publications Chalker, Lerner and Smith pre- 
sented a treatment on level distributions in disordered 
systems, which in contrast to earlier works takes care of 
a possible non-trivial structure of the eigenstates [|l0| . A 
central result of their work is a relation between the spec- 
tral form factor K(t) (i.e. the Fourier transformed of the 
two-level-correlation function R(s) ) and the ensemble 
averaged quantum return probability p(t) of a state ini- 
tially confined to a small volume. Using this result and 
from scaling theory that p(t) cx i~-° 2 / d , D 2 < d being 
the fractal (correlation) dimension |n), Chalker et al. § 
derived the critical spectral compressibility 
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where the anomalous diffusion exponent 77 is related to 
D 2 by V = d-D 2 g§. 
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In spite of the Poisson-like behavior of the number vari- 
ance £2 (N) , there still is a strong level repulsion for small 
level separations. This is obvious in the spacing distribu- 
tion function P(s), which deviates only in the tail from 
the Wigner surmise @,|,|,|l3| (Fig. |). Hence, the non- 
vanishing compressibility must be due to level density 
fluctuations on larger level separations. 

Up to now, Eq. (Q) could not be confirmed directly by 
numerical calculations. The numerical results achieved 
for the Anderson transition in d = 3 dimensions 
are compatible with a linear increase of the number vari- 
ance with N but their rather large numerical uncertain- 
ties do not allow to strictly exclude other dependencies. 
In particular, the slope of ^(N) for large N could not 
be determined precisely. Besides, the value of the mul- 
tifractal exponent D2 for the 3d-MIT is not known very 
accurately 0]. In case of the MIT in two-dimensional 
quantum Hall systems the only numerical simulation fl5| l 
presented so far shows no significant linear contribution 
to E 2 (iV). 

The difficulties in determining the compressibility are 
due to the necessity to investigate level separations that 
are large compared to the average level spacing. Since 
only the critical region (which is usually a small portion 
of the entire spectrum) can be used for the statistics, one 
has to go to very large system sizes to produce a sufficient 
number of critical levels. Unfortunately, the widths of the 
critical region achievable with present computer capaci- 
ties extends over no more than a few hundred levels, so 
that the results are strongly affected by incalculable fi- 
nite size effects. This can of course not be compensated 
by a large number of disorder realizations. Consequently, 
to study critical level statistics numerically another way 
has to be found to produce larger critical spectra. 

In this communication we show that by using the 
Chalker-Coddington network model jl6| for the Quan- 
tum Hall Effect the problems we just mentioned can be 
avoided. As pointed out by Fertig in a detailed semi- 
classical analysis leading to a description very similar 
to that in [[L6|, an energy dependent unitary operator 
U(E) associated to the network model offers an alter- 
native method for the numerical calculations of energy 
spectra and eigenstates [0. However, here we do not use 
U(E) for determining the real energy spectrum E n , but 
calculate the eigenvalues e Wl of U (E) itself, whereby the 
energy E is fixed to the critical value E c . The obtained 
quasi-spectrum loi is governed by the same statistics as 
the energy spectrum at E c and is therefore suitable for 
our purposes. By using this new method it is possible to 
improve the statistics considerably and to confirm that 
X = rj/2d to a high degree of accuracy. 

Our numerical approach is closely related to that of 
Edrei et al. Jlq] , where the concept of a network model 
has been used for calculating wave propagation through 
random media. Actually, the definition of network states 
and operator used here are in principle identical to those 



in However, the main difference are the boundary 

conditions. Whereas in Q open systems are treated, 
since they focused on transmission amplitudes, here we 
use closed systems in order to get information about the 
energy level distribution. 

The network-model for the Quantum Hall Effect |f6| is 
based on ideas developed in the early eighties |f9| for the 
description of the Anderson transition in terms of scatter- 
ing theory. It provides a semiclassical description of a 2d 
electron in a quantizing magnetic field and a smooth dis- 
order potential with correlation length A large compared 
to the magnetic length l c . The electron executes a fast 
cyclotron motion on a circle of radius l c around a guiding 
center, which drifts slowly along a contour r of constant 
energy, V{r) — E = E' — hui c /2. At saddle-points of the 
potential with energies close to E the electron tunnels 
with an appreciable probability between different con- 
tours. 

The motion of electrons along the contours is depicted 
by one-dimensional, unidirectional channels, called links. 
The electron tunneling between them is represented by 
2x2 scattering matrices S — {t m i}, which connect the 
complex current amplitudes ip k , ^1 of incoming links to 
those in outgoing links ipm,' l Pn (Fig. |]). Due to the 
random length of the links between neighboring saddle- 
points an electron acquires random phases 4>ji which we 
absorb into the scattering coefficients t m i. 

The coefficients t m i depend on the electron energy 
E and can, in principle, be determined by semiclassi- 
cal methods for a given disorder potential jl7j]. For 
a saddle-point at zero energy and tunneling energy E t 
the tunneling amplitudes are T = \t mk \ 2 = \t n i\ 2 = 
(1 + exp(-£;/£; t ))~ 1 and R = \t ml \ 2 = \t nk \ 2 = 1 — T. 
Moreover, the random phases (f>mi = arg(imz) depend in 
a rather complicated manner on the energy E |l7| . 

A state VP of the network is given by its complex am- 
plitudes on the links, W — {ijjj}j- It is stationary at 
energy E if the scattering condition 

ipm = t mk (E)iJj k + t m i(E)Tpi (3) 

is satisfied at each saddle point. 

Defining an unitary operator U(E) by 

U{E) ei =t ml {E) 

&rn ~T" tnl 

(E)e n , (4) 

where = {Sji}j, this condition can be written as 

U(E)y = (5) 

The energy enters parametrically via the coefficients 
t m i(E). This equation has non-trivial solutions only for 
discrete energies E n . According to Fertig |L7]], these 
energies E n are eigenenergies of the system and their 
eigenvectors ^> n determine the amplitudes of the corre- 
sponding eigenstates on the different equipotential con- 
tours (links). Hence, under certain circumstances the 
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scattering condition (g) offers an alternative method of 
determining numerically eigencnergics and states of an 
electron in an disordered system, as it was first — to the 
best of our knowledge — pointed out by Fertig. This 
method was utilized in p6[ for the numerical calculation 
of critical eigenstates in quantum Hall systems. (Like in 
p8[ , Eq. (^) can also be taken as a definition of a time 
evolution for states at energy E, ^(t + r) = U(E)^(t), 
whereby the energy dispersion is neglected pl"| , pO| .) 

Before we proceed it is necessary to consider the length 
and energy scales determining the critical region. The 
averaged level spacing A is related to the system size L 
and magnetic length l c via A ~ T(l c /L) 2 , T being the 
width of the disorder broadened Landau band. Close 
to the critical energy E c = the localization length £ 
is £ ~ X\E/T\- V . Therefore, the width of the critical 
energy region is A c ~ T{\/ L) 1 / 1 ' . Consequently, the 
number N c of critical levels behaves like N c = A c /A ~ 
(A/U^L/A) 2 - 1 /-. 

Note that even for a fixed ratio L/X the number N c 
can be enhanced by increasing the ratio X/l c . We em- 
phasize that in the Chalker-Coddington model this ratio 
is arbitrarily large, so that the number of critical levels is 
not restricted, which makes the model in principle very 
convenient for investigations of critical level statistics. 
However, the calculation of the energies E n by solving 
the non-linear Eq. (||) is a difficult numerical task and 
not suitable for practical purposes. 

Therefore, in order to determine spectral statistics let 
us discuss instead of Eq. (||) the eigenvalue problem 
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For a given energy E the unimodular eigenvalues e 1 ^'^ 
define quasi-energies uji(E), I = 1,...,M = dimU(E). 
They are smooth functions of the energy and do not cross 
each other, following a theorem by von Neumann and 
Wigner J23|. According to Eq. (||), the intersection points 
of the curves oji(E) with the lines ui — 0, ±2tt, ±4ir, . . . 
determine the energy levels E n . 

The flow of the levels uji(E) obeys two symmetries: 
First, the intersection points E' n with shifted lines ui' = 
O, f2±27r, . . . must exhibit the same statistics as the orig- 
inal spectrum E n , since the corresponding transformed 
operator U'{E) = e~ in U(E) belongs to the same uni- 
versality class as U(E) (W deviates from U only by a 
global phase shift, which has no influence on the statisti- 
cal properties). Second, as long as the critical regime is 
not left, \E\ < A c , the statistical properties of the flow 
uji(E) can not change significantly with energy, because 
such a change would be accompanied by a new energy 
scale inbetween A and A c E t <C T), which makes no 
physical sense. 

Due to this homogeneity in both directions and due 
to the strong repulsion of the uji(E) they must behave 
as depicted in Fig. 0: The average slope of the curves 



varies neither strongly with the level number I nor with 
the energy E. Further, this homogeneity implies that the 
intersection points ujf with a cut c crossing the band of 
curves show essentially the same statistics, independent 
of the precise position of c. For this reason, instead of 
the real energy spectrum E n one can also use a quasi- 
spectrum oji(E) with E within the critical regime for an 
analysis of the critical level statistics pi| . A big advan- 
tage of this method is the simple fact that the oji(E) are 
far better numerically accessible than the real energies 
E n . They can be calculated by solving the linear eigen- 
value problem (j|) with standard numerical methods. ]25| ] 

For our calculations we used closed networks of 50 x 50 
saddle-points with periodic boundaries in one and reflect- 
ing in the other direction. The transmission amplitudes 
were set to the critical value T — 1/2 for models describ- 
ing the transition point and to T± = (1 + cxp(±E / Et))^ 1 
with E/E t — 10 for non-critical systems with strongly 
localized states. The disorder is represented by ran- 
dom scattering phases <j) m i — arg(i m ;). Note that the 
calculations are done at constant energies E = E c and 
E = 10 E t , respectively, hence we did not have to know 
the energy dependence of the phases 4>i m (E). From this 
settings we obtained random network operators at the 
critical point, U(E = Ec), and deep in the localized 
regime, U(E = l0E t ), of dimension M = 2 x 50 x 50. Di- 
agonalizing them by standard numerical methods yields 
critical (E — E c ) and non-critical (E = 10E t ) quasi- 
spectra of M — 5000 levels each. 

As explained in the considerations given above, the 
critical level statistics can be determined by analyzing 
the critical quasi-spectra. Although the statistics do not 
change within the quasi-spectra one has to take into ac- 
count that the total number of quasi-levels per spec- 
trum is fixed to M. So, when calculating the number 
variance £2^) one has to confine oneself to interval 
sizes Aw with averaged level number N = (n)^ small 
compared to M. We checked by numerical simulations 
with Poisson distributed levels that at a total number of 
M = 5000 levels deviations from the expected number 
variance ^(N) — N are negligible for N < 300. 

For the determination of the critical number vari- 
ance E2 (N) we divided the quasi-spectra of 40 different 
disordered critical network operators U(E C ) into non- 
overlapping intervals of length Auj = (27r/5000)A, N 
ranging from 1 to 300. For each N this results in an 
ensemble of mjv = 40 x M/N intervals with level num- 
bers Hi and (rii) = N, from which we calculate the level 
number variance £2 (-A0 = 1W_ Y^h=i ( n i ~ ^) 2 - 

The results plotted in Fig. show a clearly linear be- 
havior of the number variance ^(N) for a wide range of 
N. Eq. (g) predicts a slope of r]/2d = (2 - D 2 )/2d = 
0.125 ± 0.01, where we have taken D 2 = 1.5 ± 0.05 
from independent numerical calculations of critical states 
p7| , p2 21 1 . A least square fit of our data yields a slope 
0.124 ±0.006, which agrees with the prediction very well. 
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The dashed lines mark the range of the expected fluc- 
tuations of £2 (n) due to the finite number ijin of inter- 
vals, calculated via the Xa-distribution for a = 0.8,0.2. 
This clearly indicates that the deviations of the data from 
the straight line are not systematic but due to statistical 
fluctuations. 

The spacing distribution P(s) plotted in the inset of 
Fig. |U shows for small spacings s a WD-type behavior 
for a GUE ensemble (dashed line), P(s) oc s 2 , indicating 
strong level repulsion for small level spacing. 

We use the same procedure as for the critical quasi- 
spectra for five non-critical quasi-spectra at E = 10E t 
in the strongly localized regime. Here ^{N) follows a 
straight line of slope 1, as it should be, since in this region 
the levels are Poisson-distributed. 

To summarize, the recently derived relation between 
the multifractal exponent ?y = d — D2 of cigenstates 
and the spectral compressibility at the mobility edge, 
X = T}/2d, has been confirmed numerically for the in- 
teger quantum Hall derealization transitions. This has 
been done by introducing a new method to investigate 
spectral properties of disordered systems. 
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FIG. 1. The Chalker-Coddington network. At each sad- 
dle point a scattering matrix S describes the transition from 
incoming to outgoing states. The operator U maps each in- 
coming link amplitude to the two outgoing links according to 
the transmission coefficients tmi,t„i, .... 



FIG. 2. The eigenvalues exp(iu>i(E)) of U(E) as func- 
tions of the energy E. The intersections with the lines 
Lu — 2n z,z G Z determine the energy- levels E n , those with 
the E = U c -line the quasi-energies uji{E c ). 
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FIG. 3. The level number variance E2 plotted against the 
average number of levels N at the critical point (O) and in 
the localized regime (□). The solid line has a slope of 0.124, 
the dashed lines mark the range of the expected statistical 
fluctuations (xL a — 0.8, 0.2). The dotted line has a slope of 
1. The inset shows the level spacing distribution P(s) at the 
critical point (O) and the WD-distribution for GUE (dashed 
line) . 
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